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Abstract 

Strongly-interacting systems consisting of particles that interact through a large scattering length 
satisfy universal relations that relate many of their central properties to the contact, which measures 
the number of pairs with small separations. We use the operator product expansion of quantum field 
theory to derive the universal relation for the inelastic 2-body loss rate. We verify the universal 
relation explicitly by direct calculations in the low-density limit at nonzero temperature. This 
universal relation can be tested experimentally using ultracold quantum gases of atoms in hyperfine 
states that have an inelastic spin-relaxation channel. 
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I. INTRODUCTION 



Atoms whose scattering length is large compared to the range of their interactions have 
universal properties that depend only on the large scattering length and are insensitive to 
details of their interactions at short distances jl|. These properties are universal in the sense 
that they apply to other types of nonrelativistic particles with short-range interactions, such 
as hadrons and nuclei, provided those particles have the same symmetries and mass ratios 
and large scattering lengths. The strong correlations associated with the large scattering 
length present challenges to understanding and calculating the properties of many-body 
systems of such particles. These challenges are particularly acute in the unitary limit in 
which the scattering length is infinite. This limit is characterized by scale invariance, so 
there are no well-defined quasiparticles upon which to base our understanding of the system. 

Relations between various properties of a system that depend only on the large scattering 
length are called universal relations. Such relations provide powerful constraints on the 
behavior of the system that must be respected in spite of the strong correlations. Shina Tan 
derived a number of universal relations for the strongly-interacting two-component Fermi 
gas, which consists of two types of fermions that interact through a large scattering length 
IJ]. These relations apply to any state of the system: few-body or many-body, homogeneous 
or trapped, normal or superfluid, balanced or polarized. For a few-body system, they apply 
as long as the energies relative to the appropriate threshold are small compared to the energy 
scale set by the range of the interactions. For a many-body system, they apply as long as 
the temperature and number density are small compared to the energy and density scales 
set by the range. 

Tan's universal relations all involve a property of the system C called the contact that 
measures the number of pairs of atoms of the two types with small separations. The contact 
is an extensive quantity that can be expressed as the integral over space of the contact 
density: C = J d^rC. A convenient operational definition of the contact is provided by 
Tan's adiabatic relation, which gives the rate of change of the total energy E of the system 
as the scattering length is varied jsj: 

' -E^-^C, (1) 



da ^ Stt/x 

where n is the reduced mass of the two types of atoms. If the system is at nonzero tem- 
perature, the derivative in Eq. ([T]) should be taken with the entropy held fixed. Using the 
adiabatic relation, any theoretical method that can be used to calculate the total energy of 
the system can also be used to calculate the contact. If the total energy can be measured in 
an experiment and if the scattering length can be controlled experimentally, for example by 
using a Feshbach resonance, the adiabatic relation can be used to measure the contact [5|. 

Tan derived his universal relations from the many-body Schrodinger equation by using 
generalized functions designed to implement the Bethe-Peierls boundary conditions associ- 
ated with the large scattering length 0-0]- They were rederived by Braaten and Platter 
using the machineryof quantum field theory, specifically renormalization and the operator 
product expansion [6|. Tan's universal relations make it clear that the contact plays a cen- 
tral role in the strongly-interacting Fermi gas. Subsequent work revealed that several of the 
most important experimental probes of ultracold atoms are related to the contact, includ- 
ing radio-frequency spectroscopy 0, [sl and photoassociation [1, 10|. A review of universal 



relations for fermions with a large scattering length was presented in Ref. [11 
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One of the universal relations given in Ref. is for the inelastic 2-body loss rate for 
fermions with a large scattering length. If atoms of types 1 and 2 have inelastic scattering 
channels with a large energy release, low-energy atoms will be lost from the system through 
inelastic collisions. Because of the inelastic scattering channels, the scattering length a must 
be complex with a negative imaginary part. The rate at which the number of either type of 
atom decreases is proportional to the contact:^ 



The fact that the loss rate is proportional to C was first pointed out by Tan [12|. The 
prefactor of C in Eq. ([2]) was first written down by Braaten and Platter [6]. To the best of 
our knowledge, a derivation of Eq. ([2]) has not been presented previously. In this paper, we 
derive the universal relation for the 2-body loss rate using the operator product expansion 
of quantum field theory. 



II. INELASTIC TWO-BODY LOSS 

Inelastic two-atom scattering channels with a large energy release cause the loss of low- 
energy atoms from the system. The inelastic scattering channel may consist of lower hy- 
perfine states of the atoms of types 1 and 2. The loss of low-energy atoms is actually due 
to their transformation into high-energy atoms. If the system is in a trapping potential, 
the high-energy atoms may have enough energy to escape from the trap. Even if they do 
not escape, their effects on the low-energy atoms can be neglected if their interactions with 
the low-energy atoms are weak enough. If the energy gap -Egap between the threshold for 
atoms 1 and 2 and the threshold for the final-state atoms is much larger than the low en- 
ergies of interest, the inclusive effects of the inelastic collisions can be taken into account 
through anti-Hermitian terms in the Hamiltonian. These terms must be anti-Hermitian to 
account for the loss of probability from transitions of pairs of low-energy atoms to pairs of 
high-energy atoms, which are not taken into account explicitly. 

An important consequence of the large energy gap -Egap is that there must be a small 
separation between the atoms of types 1 and 2 that disappear. By conservation of energy, 
the high-energy atoms created by the inelastic scattering of atoms of types 1 and 2 will 
emerge with large relative momentum that is approximately (/ii?gap)^^^. By the uncertainty 
principle, the localized wave packet associated with each outgoing atom can be traced back 
to a region from which the atom emerged whose size is approximately h/dxEg^pY^'^- But 
the two regions from which the two final-state atoms emerge must also be separated by a 
distance at most of order h/dxEgs.pY^'^, because the interactions that change the spin states 
of the atoms must also deliver a momentum kick of order {/jEg^pY^'^. This large momentum 
can be transferred to the two final-state atoms only if they are close together. Even if the 
momentum is transferred through an intermediate atom in the original state 1 or 2, that 
atom must be far off its energy shell and therefore can propagate only over a short distance. 



^ The relation actually stated in Ref. [ll| is that the rate of change of the number density is proportional 
to the contact density C, which is not strictly correct. Integration over space is necessary to eliminate a 
term involving the divergence of the probability current. 
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The two atoms that disappear must therefore have a separation of order (fiEgi^pY/^. Thus 
the non-Hermitian terms in the Hamiltonian must have short range. 



III. DERIVATION OF THE UNIVERSAL RELATION 

In this section, we derive the universal relation in Eq. (|2]) between the inelastic 2-atom 
loss rate and the contact. We use quantum field theory methods similar to those used in 
Ref. [6[ to give concise derivations of other universal relations. 



A. Quantum field theory 

We begin by presenting the quantum field theory formulation of the problem of atoms 
with a large real scattering length a. We introduce quantum field operators ipii^r) and 
'4)2{f) that annihilate atoms of type 1 and 2, respectively. If either atom is a boson, the 
fields satisfy the canonical commutation relations. If they are both fermions, they satisfy 
anticommutation relations. The masses of the atoms are mi and m2, and their reduced mass 
is /i = 17111712/ {mi + 777.2). The Hamiltonian operator for the system of atoms is the sum 
of a kinetic term, an interaction term, and possibly a term describing an external trapping 
potential: H = T + U + V . The kinetic and interaction terms can be expressed as 

cr 



U = ^ Jd'r^lijl^^Mr), (3b) 

where go is the bare coupling constant, which depends on the ultraviolet momentum cutoff 
A. The products of quantum field operators followed by the argument (r) indicate that the 
operators are all evaluated at the same position r. In the fermionic case, the ordering of the 
fields in the interaction term guarantees the positivity of the local operator. The ultraviolet 
cutoff can be implemented by limiting the Fourier transforms of the quantum field operators 
to wavenumbers satisfying \k\ < A. The scattering length for a pair of atoms of types 1 and 
2 will have the value a if the bare coupling constant is chosen to be 

Ait 

= l/a-2A/vr- 

We now consider the problem of atoms that have inelastic scattering channels with a large 
energy gap. The large energy gap ensures that the anti-Hermitian terms in the Hamiltonian 
density must have a short range. We can therefore take them to be local operators involving 
the quantum fields ipi and ip2- Since the two atoms that disappear must have a very small 
separation, the appropriate operator is ipiip2'ip2'ipi, which annihilates a pair of particles at a 
point and then creates them again. This is precisely the operator in the interaction term of 
the Hamiltonian in Eq. ( 13bp . To take account the loss of atoms, we can add such a term 
to the Hamiltonian with a coefficient that is pure imaginary and negative. Equivalently, we 
can simply take the bare coupling constant Qq in Eq. ( ]3bp to have a negative imaginary part. 
The inelastic scattering channel gives the scattering length a a negative imaginary part. 
The complex bare coupling constant go can be obtained simply by inserting that complex 
scattering length a into Eq. (jl]). 
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B. Particle loss rate 



We now consider the loss rate of the atoms. The time evolution of a state \X{t)) in the 
quantum field theory is given by 

\X{t)) = exp{-iHt/h)\X{0)). (5) 

Because H has a non-Hermitian part, the norm of the state \X{t)) decreases with time. The 
average number density of atoms of type a in the state \X(t)) is given by the expectation 
value of the number density operator: 

Inserting the time evolution from Eq. ([5]) and differentiating with respect to time, we get 

jnAr,t) = (^{H^ £Mr) - £Mr) H) - {£Mr)){H^ - H)) • (7) 

The contribution from the kinetic term in the Hamiltonian can be expressed as the divergence 
of a current: 

([f,^t^.(r)]>=zW- (Mr)), (8) 
where J^r is the probability current operator for atoms of type a: 

Mr) = (^lV^.(r) - V^t^.(^)) . (9) 

The flow term in Eq. ([8]) can be eliminated by integrating over all space and using the 
divergence theorem. The resulting expression for the rate of change of the total number of 
particles of type a is 

jN.= 'j_ Jd'r ((f/t^t^.(r)-^t^,(r)f/)-(^t^.(r))(m-f/>), (10) 

where If is the interaction energy operator in Eq. fl3bl) . Using the (anti) commutation rela- 
tions for the field operators, U can be expressed in terms of the number density operators 
iplipi and ip2'^2- It therefore commutes with iplip^^r). Thus only the non-Hermitian part of 
U contributes to Eq. (ITU]) . The prefactor in W — IJ can be expressed as 

Thus the loss rate in Eq. f lTU]) reduces to 



go\'44^2Mr')){i^iMr))) ■ (12) 



dt 4:71(1 



We have combined the factor of [(^op with the operator ijj\ilj\ilj2ilJi, because the dependence 
of matrix elements of this operator on the ultraviolet cutoff A is exactly cancelled by the 
A-dependence of |(7oP 0- The loss rate has been expressed in terms of a double integral 
over space of the connected part of the correlation function of two operators. 
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C. Operator product expansion 



We proceed to use the operator product expansion to reduce the expression for the loss 
rate to single integrals over space of expectation values of local operators. If the separation 
|r' — rl of the operators in the first term in the integrand in Eq. flT^ is much larger than the 
correlation length of the system, the expectation value of the operator product factors into 
the product of the expectation values and it is therefore cancelled by the second term. Thus 
the integrand goes to zero as |r' — r| oo. This motivates the expansion of the operator 
product in powers of | r' — r | : 

\goM4ij2Mr')i^iMr) = S^r' - r) \go\'44ij2Mr) + J^/anlr' - r) Or^ir). (13) 

n 

This operator equation is the operator product expansion. The sum is over local operators 
On with increasingly higher scaling dimensions. The coefficients /o-n are functions of the 
separation vector r' — r that depend on the scattering length a and the mass ratio mi/m2- 
The coefficient of the leading operator is a Dirac delta function of r' — r. It can be derived 
simply by using the (anti) commutation relations for the field operators. The coefficients of 
higher dimension operators can be calculated using diagrammatic methods [g], [l^. Inserting 
the operator product expansion into Eq. (fT2|) . the loss rate reduces to 



dt' 




(14) 



The expectation value in the first integral on the right side can be identified as the contact 
density 

C(r) = {\go\'44i^2Mr))- (15) 

The additional terms involve integrals over space of expectation values of the local operators 
On- Their coefficients Jd^r' frjn{f') are functions of a and the mass ratio. Omitting the 
contributions in Eq. (HM from higher- dimension operators, we obtain the universal relation 
for the inelastic two-atom loss rate in Eq. (jj]). 

The operator product expansion (OPE) in Eq. f|T3l) is a systematic expansion in local 
operators with increasingly higher scaling dimension. The scaling dimension is a property 
of an operator that determines the behavior of its correlation functions at short distances. 
The correlation function {0{r) 0'{r')) for a pair of local operators generally diverges as 
a power of the separation |r — r'l as the separation vector approaches 0. The power is 
determined by the scaling dimensions of the operators. The quantum field operators ip^ 
and the number density operators V'JV^o- have the naive scaling dimensions | and 3 of a 
nonrelativistic free field theory. A gradient in an operator increases its scaling dimension 
by 1. For example, the kinetic energy operators V-?/^^ ■ V^^o- have scaling dimension 5. A 



complete set of dimension-5 operators is listed in Ref . jl5| . In a free field theory, the scaling 
dimension of a composite operator is the sum of the scaling dimensions of its factors. Thus 
the naive scaling dimension of ipl'^l'^'i'^i 6. However scaling dimensions can be strongly 
modified by interactions. In the case of a large scattering length, the scaling dimension 
of is reduced from 6 to 4 14 1. Three-body operators, which contain 3 quantum 

fields and their Hermitian conjugates, have non-integer scaling dimensions 1^. In the case 
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of fermions with equal masses (mi = 7712), the lowest-dimension 3-body operators that can 
appear in the OPE in Eq. (fT3|) are Ol ■ O^, where = 2'ip2ilJi'Vijja — ijJa'^ iip2i'i) and a is 1 
or 2. They have scaling dimension 8.54 and the next higher operator has scaling dimension 
9.33 [l5| . In the case of unequal masses mi < m2, the operator analogous to 0{ ■ Oi has 
a scaling dimension that decreases as 1112/1711 increases. It reaches 5 at the critical value 
1712/1711 = 13.6 and remains at 5 for larger values of m2/mi 1^. If atom 1 is a boson, the 
lowest- dimension 3-body operator that can appear in the OPE in Eq. ([13]) is ?/'iViVI^2V^i^i. 
The scaling dimensions of this operator is 5 for all values of the mass ratio 1112 /mi. 

We now consider which higher- dimension operators appear in the OPE for the loss rate 
in Eq. (fT^ . The only operators that can contribute are those for which the integral over 
all space of the expectation value ((9„(r)) is nonzero. Operators that are total derivatives 
have expectation values whose integrals are zero. If the system is homogeneous, only scalar 
operators have nonzero expectation values. If the system is static, operators that are total 
time derivatives have expectation values that are zero. Linear combinations of operators 
that vanish upon using the equations of motion for the quantum field operators also have 
zero expectation values. The additional operators that appear in the OPE in Eq. f|T^ are 
also limited to 3-body operators. This can be seen by using the (anti) commutation relations 
to normal order the quantum fields in the operator product in Eq. f|T3l) . The result is the 
sum of the first term on the right side and the normal-ordered form of the left side. But 
the normal-ordered operator has nonzero expectation value only in states that include 3 or 
more particles. Thus only 3-body and higher-body operators can appear in the additional 
terms in the loss rate in Eq. ([T^ . The contact operator in Eq. f[T^ has scaling dimension 4. 
In the case of equal-mass fermions, the next term has an operator with scaling dimension 
8.54. In the case of fermions with m2/mi > 13.6 and in the case of bosons, the next term 
has an operator with scaling dimension 5. 

We next determine how the coefficients Jd^r' fcmir') of the higher- dimension terms in the 
loss rate in Eq. f|T^ depend on a. The normalizations of the operators On can be chosen so 
that their matrix elements have nontrivial unitary limits as a — )■ ±00. A good example is 
the contact density operator |(yfoPV'iV'2'V^2'?/'i; for which the prefactor l^fQp guarantees a well- 
behaved unitary limit. Given such operators, we can use dimensional analysis to determine 
how their coefficients depend on a. The first term inside the parentheses in Eq. (1141) is 
the contact, which has dimensions 1/Length. If the operator On has scaling dimension 
A„, jd?r{On{r)) has dimensions (Length)^"^". By dimensional analysis, its coefficient 
JcPr' fcjn{f') niust have dimensions (Length)^""^. Since a is the only dimensionful variable 
the coefficient depends on, Jd^r' f^^ni''^') must be |a|^"~^ multiplied by a numerical constant 
that depends only on the sign of a and on the mass ratio 7712/ mi. 

The truncation of the OPE for the loss rate in Eq. (|T^ to the leading term as in Eq. ([2]) is 
obviously a good approximation in the weak-coupling limit a — )■ 0. Let the largest important 
momentum scale for the system be hk^. Its contact is then of order /cq. If On has scaling 
dimension A„, Jd^r{On{r)) will be at most of order k^"~^. Multiplying by its coefficient, 
which is proportional to |a|^"~^, we see that its contribution to the loss rate is suppressed 
by at least a factor of (A;o|a|)'^"~'^. Thus, in the weak-coupling limit, the contributions from 
higher dimensions operators are parametrically suppressed by powers of /col*^!- In the case 
of equal-mass fermions, the contribution to the loss rate from the first higher-dimension 
operator is suppressed by a factor of (A;o|a|)'^'^'^. 

There are good reasons to believe that the truncation of the OPE for the loss rate in 
Eq. (IT^ remains a good approximation even near the unitary limit a — )■ ±00. There is no 
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parametric suppression of higher dimension operators in this hmit, so the vahdity of the 
truncation depends on the convergence properties of the OPE. The OPE was discovered 



independently by Kadanoff, Polyakov, and Wilson in 1969 [16Hl8[. Wilson proposed that 
the OPE was an asymptotic expansion in the small separation |r — r'| of the operators [l8| . 
Studies of the convergence of the OPE in relativistic quantum field theories suggest that it 
is actually a convergent expansion. Mack showed in 1977 that the OPE for a conformally- 
invariant field theory in 4 space-time dimensions converges for finite | r — r' | , at least when 



acting on the vacuum state [19|. More recently, Pappadopulo, Rychkov, Espin, and Rattazzi 
showed that the OPE for a conformal field theory in any space-time dimension not only has 
infinite radius of convergence in the separation | r — r' | , but the convergence at fixed | r — r' | 
is actually exponential in the scaling dimension j20[. If the OPE is truncated by omitting 
terms with scaling dimension Aq and higher, the error decreases like exp(— AAq), where A 
depends on the separation. Furthermore, Hollands and Kopper showed that the OPE for the 
simplest interacting scalar field theory in 4 space-time dimensions has similar convergence 
properties order by order in perturbation theory [21]. It has infinite radius of convergence 
in the separation |r — r'|, and the convergence at fixed |r — r'| is exponential in the scaling 
dimension. If the OPE is truncated by omitting terms with scaling dimension Aq and higher, 
the error decreases like A^"/(Ao!)'^''^, where A depends on the state and is proportional to 



the separation [21 



The convergence of the OPE in relativistic field theories does not guarantee its conver- 
gence in nonrelativistic field theories. However it is at least plausible that the OPE in a 
nonrelativistic field theory has infinite radius of convergence in the separation | r — r' | , and 
that the convergence at fixed |r — r'| is exponential in the scaling dimension. This suggests 
that the OPE for the loss rate in Eq. (1141) may have exponential convergence in the scaling 
dimension at any fixed scattering length, including the unitary limit a = iboo. In the case of 
a balanced gas of equal- mass fermions at zero temperature near the unitary limit, the only 
important momentum scale is hkp, where kp = (Svr^n)^/^. The contribution to the loss rate 
in Eq. f lT^ from an operator of scaling dimension A^ is not parametrically suppressed, but 
instead is enhanced by a factor of (/c^lal)^""^. However if the convergence is exponentially 
fast in the scaling dimension, the exponential factor exp(— AA^) overcomes the parametric 
enhancement. Thus the universal relation in Eq. obtained by truncating the OPE after 
the leading term may a good approximation even in the unitary limit. 



IV. CONTACT IN THERMAL EQUILIBRIUM 

In this section, we calculate the contact for a single pair of atoms in thermal equilibrium. 
We then use the result to obtain the contact density for a dilute gas of atoms or dimers. 
In this section, we take the scattering length a to be real and ignore the effects of inelastic 
scattering channels. 



A. Single pair of atoms 

We consider two atoms that are restricted to a volume V that is much larger than \a\^ and 
are in thermal equilibrium at temperature T. The adiabatic relation in Eq. ([T]) can be used 
as an operational definition of the contact. It can be generalized to nonzero temperature by 
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taking the derivative with the entropy fixed. Equivalently, it can be expressed in the form 



(16) 



da~^ Svr/i 

where F is the free energy and the derivative is taken with the temperature fixed. The free 
energy can be determined from the partition function: 

Z = exp(-/3F), (17) 

where /3 = 1/(/cbT). We will use the universal relation in Eq. (fT6l) to calculate the contact 
for the single pair of atoms. 

In the center-of-mass frame, the spectrum of the pair includes a continuum above the 
threshold from scattering states. Their energies can be parameterized as E = /i^/c^/(2//), 
where k is the wavenumber. The phase shift for an S-wave scattering state is 

5{k) = — arctan(A;a). (18) 

If a > 0, the spectrum also has a discrete energy below the threshold: the dimer, whose 
energy is = —Ed, where E^ is the binding energy: 

Since the center-of-mass motion and the relative motion of the pair decouple, the partition 
function is the product of the associated partition functions Zcm and Zj.e\. The dependence 
on a enters only through Z^-ei. In the limit a — >■ 0~, there are no interactions and Zrei reduces 
to the partition function for a free particle. 

where V is the spatial volume and At is the thermal wavelength: 



We express Zy.^\ as the sum of the free partition function in Eq. (!20|) and an interaction 
term that has all the dependence on a: Zroi = ^froe + Z\^f The interaction term receives 
contributions from the continuum of S-wave scattering states and also, if a > 0, from the 
dimer. The spectrum of scattering states can be made discrete with energies E^ by making 
the spatial volume compact. A convenient method for discretizing the scattering states is 
to put the system in a spherical box with a large radius R. The interaction term can then 
be expressed as 



^int(a) = ^(a)e^^- + V g-'^^'^-/^^^) - V e-f^^'^H^''^) . (22) 

n=l n=l 

The scattering wavefunctions sin[fcr + 6{k)] must satisfy the boundary condition 

sm[knR + 6{kn)] = 0. (23) 
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Its nontrivial solutions satisfy 



knR — arctan(/c„a) = nvr, n = 1,2, 3, ... . (24) 

In the limit i? — )■ oo, the discrete energies approach a continuum and the sum over integers 
approaches an integral: 

oo 



TT 

n=l 



^ 1 r I \ 



Inserting this into Eq. fl22|) . we find that the dependence on R cancels and the interaction 
term reduces to 

= ^(a)e'^'^'/(2^'^') - - Tdk ^—^e-^^''''/^'^\ (26) 

TT Jq 1 + fc O 

The universal relation in Eq. ( !T6|) implies that the contact for the single pair of atoms is 
given by 

Cpair = ^^-1 l°S(^free + ^int)- (27) 

The derivative of Z^^t is evaluated most easily if we first change the integration variable in 
Eq. fl26|) to the dimensionless variable x = ak. The resulting expression is 



^""^ Z,^, = e{a)-e^^'/^^^'^'^ + ^a' H dk ^^e-^^''''/^'^\ (28) 



h^(3da-^ ' ' a Jo l + a^F 

Our final result for the contact for a single pair of atoms is 

_ 8nfi{d/da-^)Zi^t /^qv 

where Zjnt and its derivative are given in Eqs. ( 126|) and fl28|) . The integrals in these expres- 
sions can be evaluated analytically in terms of error functions. 

The contact in Eq. (1291) is for a pair of atoms in thermal equilibrium and also in chemical 
equilibrium with respect to transitions between a dimer and two unbound atoms. The 
contact for two unbound atoms that are not in chemical equilibrium with the dimer can be 
obtained by omitting the dimer contributions in the numerator and denominator of Eq. fl29|) . 
In the large volume limit, we need only keep the first term in the denominator. The contact 
for the two unbound atoms reduces to 

The contact for a dimer that is not in chemical equilibrium with unbound atoms can be 
obtained by keeping only the dimer terms in the numerator and denominator of Eq. (!29|l . 
The result is extremely simple: 

Co = (31) 
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B. Dilute gas of atoms or dimers 



We can use the expression in Eq. (1301) for the contact for a single pair of unbound atoms 
to deduce the contact density for a homogeneous gas of atoms in the dilute limit. If the 
numbers of atoms in the volume V are A^^i and N2, the total number of pairs is N1N2. Each 
pair gives the same contribution Caa to the contact. We obtain the contact density by 
dividing by the volume: Caa = NiN2Caa/V- Using the expression for Caa in Eq. f l5U]) . we 
find that the contact density is proportional to the product of the number densities ni and 
^2: 

Caa = 8a^A^ (^j^^dk Y^2^~'''^'^'^"'^) ^i^^. (32) 



We can use the expression in Eq. ( 131]) for the contact Cd of a dimer to deduce the 
contact density for a homogeneous gas of dimers in the dilute limit. The total contact is 
multiplied by the number iV^) of dimers. We obtain the contact density by dividing by the 
volume. The result is proportional to the number density ud oi the dimers: 

Cd = —nD- (33) 
a 



V. DIRECT CALCULATIONS OF INELASTIC LOSS RATES 

In this section, we carry out direct calculations of inelastic loss rates for dilute homo- 
geneous gases of atoms or dimers in thermal equilibrium. We take the scattering length a 
to be complex, with a negative imaginary part that takes into account inelastic scattering 
channels. For the direct calculations, it is more convenient to use the inverse scattering 
length 7 = 1/a, which has a positive imaginary part. 



A. Low-density gas of atoms 

The elastic scattering amplitude for atoms with a large complex scattering length = 1/7 
at wavenumber k is 

f{k) = -^—T . (34) 
—7 — ik 

The differential cross section \ f{k)\^ is isotropic. The elastic cross section is obtained by 
multiplying it by the solid angle: 



cr 



(elastic) 



47r 



(^) = 1 • (35) 

|— 7 — ik\ 



According to the optical theorem, the total cross section is the imaginary part of the forward 
scattering amplitude multiplied by An/k: 



k \—1 — if 



(36) 
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The inelastic cross section is the difference between the total and elastic cross sections in 
Eqs. (l36l) and (135|) . The inelastic transition rate g{k) per volume squared can be obtained 
by dividing the inelastic cross section by the ffux factor fi/hk: 

9W = — (37) 

In the dilute limit, the inelastic loss rate can be written 

d d , , 

— m = — n2 = -K2nin2 , (38) 
dt dt 

where rii and n2 denote the number densities of atoms 1 and 2, respectively. For an ensemble 
of particles 1 and 2, the inelastic loss rate coefficient K2 can be calculated via a statistical 
average of the inelastic transition rate g{k). Taking the ensemble to be atoms in thermal 
equilibrium at temperature T, the rate coefficient is 

Inserting the inelastic transition rate g{k) in Eq. f l37|) . we obtain the loss rate coefficient K2'- 
K2 = ^ Im(7) rffc ^^e-^^'''/^'^'^ . (40) 



According to the universal relation in Eq. ([2]), the right side of Eq. ( I38l) must be propor- 
tional to the contact density: 

K2nin2 = Im(l/a) Caa- (41) 
47r/i 

Comparing the expression for K2 in Eq. ( HOl) and Caa in Eq. ( I32l) . we see that they agree 
only to leading order in Im(a). This is not surprising, because the contact density in Eq. f l32ll 
was calculated only for the case of real a. That calculation could presumably be extended 
to the case of complex a. The result must agree with Eq. fHTl) in the dilute limit. Thus the 
contact density in the dilute limit for atoms with a complex scattering length a must be 

Caa = 8A^ I / dk ^_y^_^j^^2 I ^1^2 • (42) 



B. Low-density gas of dimers 

The inelastic loss rate for a dilute, noninteracting gas of dimers can be written as 

d 

—ud = -Tdtid, (43) 

where F^) is the decay rate or inverse lifetime of the dimer. The decay rate can be expressed 
as an integral over momenta of the dimer wave function tp{k) in momentum space, its 
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complex conjugate il!*{k'), and the imaginary part of the two-particle-irreducible transition 
amplitude U{k, k') for a pair of atoms: 



r. = ^y j^^m{-i^u{k,k'))r{n m 

In the field theory formulation of the problem with the interaction term in Eq. f l3b[) . U{k, k') 
is simply the bare coupling constant goh!^ /{2fj). Since this is momentum independent, Tjy 
simplifies to 

n 



D = - y-lmgo) 



(A) j3l 

^ m 



The normalized dimer wave function for complex scattering length a = I/7 is 



(45) 



If the integral in Eq. ( I45p was convergent, it would be the wave function at the origin in 
position space. However it is ultraviolet divergent, and it has therefore been regularized 
with a momentum cutoff A. The regularized integral is 

^ k^dk JSTiReh) , . 

This dependence on A is precisely what is required to cancel the A-dependence of the factor 
of Im^fQ in Eq. fHSl) . The imaginary part of the coupling constant in Eq. (jl]) can be 
expressed as 

Im,o = -^^^^^. (48) 
Inserting Eqs. f HS]) and fH7|) into Eq. f HS]) . we obtain a simple result for the decay rate: 

271^ / X. / X 2;iRe(a)Im(l/a) , 

Yb = — Re(7) Im(7) = ■ (^9) 

fj, fi\a\^ 

This result for Tr, can also be obtained from the expression for the dimer binding energy in 
Eq. fll9l) by analytically continuing the scattering length a to complex values. By expressing 
the complex energy —h^/{2jj,a?) in the form —Ejj — ihT£)/2, we can read off the width HTd 
of the resonance. This prescription gives the same result for Td as Eq. fl49l) . 

According to the universal relation in Eq. ([2]), the right side of Eq. ( H3|) must be propor- 
tional to the contact density: 

TDnD = -^lm{l/a)CD. (50) 
47r/i 

Comparing the expression for Td in Eq. (149|) and Cd in Eq. (!33|) . we see that they agree 
only to next-to-leading order in Im(a). This is not surprising, because the contact density 
in Eq. ( 133|) was calculated only for the case of real a. That calculation could presumably 
be extended to the case of complex a. The result must agree with Eq. (!50|) in the dilute 
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limit. Thus the contact density for a dilute gas of dimers in the case of a complex scattering 
length a must be 

87rRe(a) 

Cd = — — riD- (51) 



a 



Kohler, Tiesinga, and Julienne have derived a universal relation between the inelastic 
rate coefficient K2 and the dimer lifetime I/Td (22|:^ 

K2 = 27ra^TD. (52) 

They considered the case in which Im(a) is negligible compared to Re(a). Using the expres- 
sions for K2 in Eq. ( HOj) and Tjj in Eq. ( |49|) . we can see that the relation holds for K2 in the 
zero-temperature limit and we can generalize it to a complex scattering length a: 

This differs from the relation for real a in Eq. (l52l) at next-to-next-to- leading order in Im(a). 

VI. TRAPPED GAS OF EQUAL-MASS FERMIONS 

The universal relation for the inelastic two-body loss rate in Eq. ( fT4l) can be tested ex- 
perimentally using atoms in hyperfine states that have an inelastic spin relaxation channel. 
Such an experiment has been carried out for identical bosons by Thompson, Hodby, and 
Wieman using an ultracold trapped thermal gas of ^^Rb atoms near a Feshbach resonance 



23( 1 . The results were interpreted theoretically by Kohler, Tiesinga, and Julienne in terms 
of the spontaneous dissociation of dimers [22|]. The calculated lifetime of the dimers agreed 
well with the results extracted from the experiment except very close to the Feshbach reso- 
nance, where the predicted lifetime was too large. The extension of our universal relation for 
the atom loss rate to the case of identical bosons should allow a quantitative understanding 
of the loss rate much closer to the Feshbach resonance. Identical bosons introduce addi- 
tional complications. One theoretical complication is that the operator product expansion 
analogous to Eq. f|T3l) includes the 3-body contact operator, which has scaling dimension 5 
[23 |. The contribution from the 3-body contact may therefore not be strongly suppressed 
compared to the leading term from the 2-body contact. An experimental complication is 
that the contribution to the loss rate from 3-body recombination into the shallow dimer 



may also be important [23|. It is possible that part of the loss rate attributed to 3-body 
recombination in the experiment in Ref. j23| was actually due to the contact from unbound 
atoms. 

An analogous experiment with equal-mass fermions would be particularly easy to inter- 
pret. The severe suppression of higher dimension operators implies that the simple universal 
relation in Eq. ([2]) should be very accurate. Additional contributions to the loss rate from 
3-body recombination into the shallow dimer are also suppressed. The experiment could 



be carried out in the same way as the experiment with Rb atoms in Ref. [23|. It would 



begin with a large number of atoms in thermal equilibrium at a scattering length on the 



^ This is the relation for fermions with two spin states. For identical bosons, the proportionality factor is 
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negative side of a Feshbach resonance. The scattering length would then be ramped to a 
large positive value. The number of dimers that are produced could be controlled by varying 
the ramp rate. The scattering length would be held at the large positive value for a variable 
time, during which atoms and dimers could disappear. Finally, the scattering length would 
be ramped back to the negative side of the Feshbach resonance to dissociate the dimers, and 
the total number of remaining atoms would then be measured. By analyzing the dependence 
of the number of remaining atoms on the holding time, one could determine the contact for 
an unbound atom pair and/or the contact for a dimer. These contacts could be measured 
as functions of the scattering length and compared with the universal predictions. 

We therefore consider a dilute gas of fermionic atoms and dimers that are trapped in a 
harmonic potential and in thermal equilibrium at a very low temperature T. We take the 
two types of atoms to be different spin states, so they have equal masses mi = m2 = m, 
and we assume that they both feel the same trapping potential: 

V{r) = ^m{colx' + coy + coy). (54) 

The number densities ni and n2 for the two spin states are Gaussians whose integrals over 
space are the total numbers of atoms A^^i and N2, respectively. The dimers feel the potential 
2V{r). Their number density ra^ is also a Gaussian whose integral over space is the total 
number of dimers Nd. 

For simplicity, we assume that the imaginary part of the scattering length is extremely 
small compared to its real part. It therefore enters only in the prefactor of the contact in 
the universal relation in Eq. f|T^ . The contact has contributions from unbound atoms and 
from dimers. The contribution from unbound atoms can be obtained by using the result for 
the contact density in Eq. fl32l) together with the local density approximation. It is obtained 
by integrating the contact density in Eq. fl32|) over all space. The result can be expressed as 
N1N2CAA, where the contact for a single pair of atoms is 

and u = (uxUiyUzY^^ is the geometric mean of the trapping frequencies. The contribution 
to the contact from dimers can be obtained by using the result for the contact density in 
Eq. f l33l) together with the local density approximation. Integrating the contact density over 
all space, we obtain the simple result NdCd, where Co = Svra/m is the contact for a single 
dimer. 

The number of dimers evolves independently from the number of unbound atoms. The 
time dependence for the number of dimers is 

Nnit) = ND{0)exp{-TDt), (56) 

where = 4^ Im (l/a)/(ma). If the system is balanced, with equal numbers of atoms in 
the two spin states, the time dependence for the number of unbound atoms is also simple: 



^ Mm(l/a) ^^^^- 



_N^{0) 27rm 



(57) 



where Caa is the contact for a pair of unbound atoms in Eq. fl55|) . If the initial number 
of dimers is sufficiently large, the loss rate will at first be dominated by the dimer contact. 
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However, as the dimers decay away, the loss rate will eventually be dominated by the contact 
from unbound atoms. The cross-over point at which the dominant loss mechanism changes 
from the dimer contact to the contact from unbound atoms is when the numbers of atoms 
and dimers satisfy N1N2/ND ~ Cd/Caa- 

VII. SUMMARY 

In summary, we have developed a rigorous theoretical framework for calculating the in- 
elastic two-body loss rate for ultracold atoms with a large scattering length. We used the 
operator product expansion (OPE) to derive the universal relation for the inelastic two-body 
loss rate. The general result is given in Eq. (I14p as an expansion in terms of integrals of 
expectation values of local operators with increasingly higher scaling dimensions. If we trun- 
cate the OPE to the leading term, we get the simple universal relation in Eq. ([2]) between 
the loss rate and the contact. In the weak-coupling limit in which the largest important 
momentum scale /cq of the system satisfies fcol^l "C 1, the OPE provides a systematically im- 
provable approximation to the loss rate, because higher dimension terms are parametrically 
suppressed by powers of /co|a|. The truncation of the OPE to the leading term is particularly 
accurate in the case of fermions with equal masses, because the first correction is suppressed 
by a factor of (/co|fl|)^'^^- If the OPE has the same convergence properties as in relativistic 
field theories, the general universal relation for the loss rate in Eq. f[T^ may remain useful 
near the unitary limit. The parametric enhancement of higher dimension terms by powers 
of /co|o| could be overcome by numerical convergence factors that are exponential in the 
scaling dimension. This could justify the truncation of the OPE to the leading term, so the 
simple universal relation in Eq. (|2]) between the loss rate and the contact may be a good 
approximation even in the unitary limit. 

We verified the universal relation by direct calculations for homogeneous gases of atoms 
and dimers in thermal equilibrium in the dilute limit. The contact density for a dilute gas 
of atoms is given in Eq. fl32l) for the case of a real scattering length and in Eq. fl42p for the 
case of a complex scattering length. The corresponding results for a dilute gas of dimers are 
given in Eqs. ( 15^ and 

We described how the universal relation in Eq. ([2]) could be tested using dilute ultracold 
gases of fermions with two spin states. The time dependence of the loss rate at a fixed 
scattering length could be used to measure both the contact for a single dimer and the 
contact for a pair of unbound atoms. The dependence of these contacts on the scattering 
length could then be compared with universal predictions. Agreement with the predictions 
would provide a beautiful illustration of how the universal relations constrain the interplay 
between few-body and many-body physics. 
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